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Abstract 

We assume that an individual invests in a financial market with one riskless and one risky asset, with the latter's price 
following a diffusion with stochastic volatility. Given the rate of consumption, we find the optimal investment strategy 
for the individual who wishes to minimize the probability of going bankrupt. To solve this minimization problem, we 
use techniques from stochastic optimal control. 
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1. Introduction 

Pension actuaries traditionally have computed the liabilities for defined benefit (DB) pension plans; however, 
more and more employees are participating in defined contribution (DC) plans. Indeed, in June 2007, the Employee 
Benefits Research Institute (EBRI) reported that in 1979, among active workers participating in retirement plans, the 
percentages in DB plans only, DC plans only, and both DB and DC plans were 62%, 16%, and 22%, respectively. The 
corresponding percentages in 2005 were 10%, 63%, and 27%, respectively. 

In terms of numbers of employees, EBRI reported that in 1980, 30.1 million active workers participated in DB 
plans, while 18.9 million workers participated in DC plans. The corresponding numbers in 2004 were 20.6 and 52.2 
active workers, respectively. Finally, in terms of numbers of plans in the private sector, in 1980, there were 148 
thousand DB plans and 341 thousand DC plans; the corresponding numbers in 2004 were 47 and 653 thousand plans, 
respectively. 

Therefore, however one measures the change in employee coverage under DB versus DC plans, it is clear that 
pension actuaries will need to adapt to the migration from DB to DC plans. One way that they can adapt is to switch 
from advising employers about their DB liabilities to providing investment advice for retirees and employees in DC 
plans. The purpose of our proposed research is to help train actuaries for this opportunity under the easy-to-explain 
goal of an employee or retiree avoiding bankruptcy. 

Previous work focused on finding the optimal investment strategy to minimize the probability of bankruptcy under 
a variety of situations: (1) allowing the individual to invest in a standard Black-Scholes financial market with a rate 
of consumption given by some function of wealth, ( Young) |2004| [Bayraktar and Young] |2007a] l; (2) incorporating 
immediate and deferred annuities in the financial market, (Milevs ky et al.| |2006| |Bayraktar and Young| |2009[) ; (3) 
limiting borrowing or requiring that borrowing occur at a higher rate than lending, ( |Bayraktar and Young 2007b) ; (4) 



modeling consumption as an increasing function of wealth or as a random process correlated with the price process 



of the stock, (IBayraktar and Young 2008 ; Bayraktar et al. , 2008 Bayraktar and Young||201 1). Throughout this body 



of work, the price process of the stock is modeled as a geometric Brownian motion, which is arguably unrealistic, but 
has given results that one can consider to be "first approximations." Here we extend some of the previous work and 
allow the stock price to exhibit stochastic volatility. Additionally, we intend to find easy-to-implement rules that will 
result in nearly minimal probabilities of bankruptcy under stochastic volatility. 
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The rest of the paper is organized as follows. In Section 2, we introduce the financial market and define the problem 
of minimizing the probability of lifetime ruin. In Section 3, we present a related optimal controller-stopper problem, 
and show that the solution of that problem is the Legendre dual of the minimum probability of lifetime ruin. By solving 
the optimal controller-stopper problem, we effectively solve the problem of minimizing the probability of lifetime ruin. 
Relying on the results in Section 3, we find an asymptotic approximation of the minimum probability of ruin and the 
optimal strategy in Section 4. On the other hand, in Section 5, relying on the Markov Chain Approximation Method, 
we construct a numerical algorithm that solves the original optimal control problem numerically. In Section 6, we 
present some numerical experiments. 

We learn that the optimal investment strategy in the presence of stochastic volatility is not necessarily to invest 
less in the risky asset than when volatility is fixed. We also observe that the minimal probability of ruin can be 
almost attained by the asymptotic approximation described in Section 5.1. Also, if an individual uses the investment 
prescribed by the optimal investment strategy for the constant volatility environment while updating the volatility 
variable in this formula according to her observations, it turns out she can almost achieve the minimum probability of 
ruin in a stochastic volatility environment. 



2. The Financial Market and the Probability of Lifetime Ruin 

In this section, we present the financial ingredients that make up the individual's wealth, namely, consumption, a 
riskless asset, and a risky asset. We, then, define the minimum probability of lifetime ruin. 

We assume that the individual invests in a riskless asset whose price at time f, X t , follows the process dX t = 
rX,dt, Xo — x > 0, for some fixed rate of interest r > 0. Also, the individual invests in a risky asset whose price at 
time t, S t , follows a diffusion given by 

dS t = S t {fxdt + o- t dB < t 1) ), S = S>0, (2.1) 

in which fi > r and <x, is the (random) volatility of the price process at time t. Here, B m is a standard Brownian motion 
with respect to a filtered probability space (£2, T , P, F = { < F t }t>o)- We assume that the stochastic volatility is given by 

cr t = f(Y t ,Z t \ (2.2) 

in which / is a smooth positive function that is bounded and bounded away from zero, and Y and Z are two diffusions. 
Below, we follow Fouque et al. ( 2003) 1 in specifying the dynamics of Y and Z. Note that if / is constant, then S 



follows geometric Brownian motion, and that case is considered by Young (2004). 

The first diffusion Y is a fast mean-reverting Gaussian Ornstein-Uhlenbeck process. Denote by 1 je the rate of 
mean reversion of this process, with < e <K 1 corresponding to the time scale of the process. Y is an ergodic 
process, and we assume that its invariant distribution is independent of e. In particular, the invariant distribution is 
normal with mean m and variance v 2 . The resulting dynamics of Y are given by 

1 Pi 
AY, = - (m - Y,) dt + vJ- ABf\ Y = y e R, (2.3) 

£ V £ 

in which B (2) is a standard Brownian motion on (£2, T, P, F). Suppose B (l) and B (T) are correlated with (constant) 
coefficient p\2 e (— 1, 1). 

Under its invariant distribution N(m, v 2 ), the autocorrelation of Y is given by 

E [(Y t - m){Y s - m)] = v 2 e -|M/e . (2.4) 

Therefore, the process decorrelates exponentially fast on the time scale e; thus, we refer to Y as the fast volatility 
factor. 

The second factor Z driving the volatility of the risky asset's price process is a slowly varying diffusion process. 
We obtain this diffusion by applying the time change f — > S ■ t to a given diffusion process: 

dZ, = g(Z,) dt + h(Z,) dB„ (2.5) 
2 



in which < 6 <k 1 a nd B is a standard Brownian motion. The coefficients g and h are smooth and at most linearly 



growing at infinity, so (2.5 1 has a unique strong solution. Under the time change t — > 8 ■ t, define Z t = Zg. t . Then, the 
dynamics of Z are given by 

dZ, = 6g(Z t ) dt + h(Z t ) dB & .„ Z = z e R. (2.6) 
In distribution, we can write these dynamics as 

dZ f = Sg(Z t ) dt + ^5h(Z,) dB^, Z = z e R, (2.7) 

in which B (3) is a standard Brownian motion on (£2, T, P, F). Suppose and B (3) are correlated with (constant) 
coefficient p^ e [—1,1]. Similarly, suppose B (T) and Z? (3) are correlated with (constant) coefficient P23 e [ _ L 1]- 
To ensure that the covariance matrix of the Brownian motions is positive semi-definite, we impose the following 
condition on the p's: 

1 + 2pi 2 pi3P23 - P? 2 - Pl3 - Pl3 ^ °- ( 2 - 8 ) 

Let W, be the wealth at time t of the individual, and let n, be the amount that the decision maker invests in the risky 
asset at that time. It follows that the amount invested in the riskless asset is W, - n,. We assume that the individual 
consumes at a constant rate c > 0. Therefore, the wealth process follows 

dW, = [rW, + r)n, -c]dt + f(Y„Z,) n, dB^\ (2.9) 

and we suppose that initial wealth is non-negative; that is, Wo = w > 0. 

By lifetime ruin, we mean that the individual's wealth reaches zero before she dies. Define the corresponding 
hitting time by to := inf{f > : W, < 0). Let tj denote the random time of death of the individual, which is 
independent of the Brownian motions. We assume that Td is exponentially distributed with parameter A (that is, with 
expected time of death equal to I /A); this parameter is also known as the hazard rate, or, force of mortality. 

|Moore and Young] ( |2006| l minimize the probability of lifetime ruin with varying hazard rate and show that by 
updating the hazard rate each year and treating it as a constant, the individual can quite closely obtain the minimal 
probability of ruin when the true hazard rate is Gompertz. Specifically, at the beginning of each year, set A equal to 
the inverse of the individual's life expectancy at that time. Compute the corresponding optimal investment strategy 
as given below, and apply that strategy for the year. According to the work of Moore and Young (2006), this scheme 



results in a probability of ruin close to the minimum probability of ruin. Therefore, there is no significant loss of 
generality to assume that the hazard rate is constant and revise its estimate each year. 

Denote the minimum probability of lifetime ruin by if/(w, y, z), in which the arguments w, y, and z indicate that one 
conditions on the individual possessing wealth w at the current time with the two factors Y and Z taking the values y 
and z, respectively, then. Thus, if/ is the minimum probability that to < tj, in which one minimizes with respect to 
admissible investment strategies n. A strategy n is admissible if it is ^-progressively measurable, and if it satisfies 
the integrability condition J Q n 1 ^ ds < oo almost surely for all t > 0. Thus, ip is formally defined by 

if,(w,y,z) = mfP K -y- z [T <T d ]. (2.10) 

7T 

Here, P lv '- V ' z indicates the probability conditional on Wo = w, Yq — y, and Zq — z- Note that if w > c/r, then if/(w,y, z) = 
because the individual can invest c/r of her wealth in the riskless asset and generate a rate of income equal to c, 
which exactly covers her consumption. Therefore, we effectively only need to determine the minimum probability of 
lifetime ruin and corresponding optimal investment strategy on the domain D := {(w,y, z) e R 3 : w e [Q,c/r]}. 



3. Computing the Minimum Probability of Lifetime Ruin 

3.1. A Related Optimal Controller-Stopper Problem 

In this section, we present an optimal controller-stopper problem whose solution is the Legendre dual of the 
minimum probability of ruin ifr. It is not clear a priori that the value function if/ is convex or smooth due to the implicit 
dependent on the initial values of the state variable. By passing to the controller-stopper problem, however, we can 
obtain the regularity of 0- in more simply, which, in turn, provides an intermediate tool in the proof of regularity of ifi. 
The dual relationship and the analysis of the controller-stopper problem are, therefore, crucial and worth investigating. 
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First, note that we can represent the three Brownian motions from Section 2 as follows: given Z? (2) , and B (3) , 
define B m , B {2 \ and B (3) via the following invertible system of equations: 



d(2) _ - r(D , F, I2~r(2) 
B, - PnB, + yi -p l2 B, , 



B? ) =p l3 BY> + 



, P23 -P12P13 r(2) 



J(l -Pi 2 )(l -P[ 3 ) - (P23 -P12P13) 2 __, 



(3.1) 



V 1 ~^2 V 1 _P 12 

One can show that Z? (1) , fi (2) , and fi (3) thus defined are three independent standard Brownian motions on (O, T , P,F). 
Also notice that condition (2.8 1 on the p's guarantees that the expression under the square root in the coefficient of 
B^ is non-negative. 



Next, define the controlled process X 7 by 



ax; = -(r-A)x;dt- 



M,z,) 



r XjdB^+y^dB^+yfdB^, 



X = x>0, 



(3.2) 



in which y = (y <2) , y (3) ) is the control, and Y and Z are given in (2.3 1 and (2.7 1, respectively. 
For x > 0, define the function $r by 



inf sup E v 

T y 



J e- A 'cX] dt + e- AT min ((c/r)X?, l) 



(3.3) 



is the value function for an optimal controller-stopper problem. Indeed, the controller choo ses among processes 
y in order to maximize the discounted running "penalty" to the stopper given by cXj in (3.3 1. On the other hand, 
the stopper chooses the time to stop the game in order to minimize the penalty but has to incur the terminal cost of 
min ( {clr)Xj, l), discount ed by e~ iT when she stops. 

Bayraktar and Young (20 11) consider a controller-stopper problem that is mathematically similar to the one in this 



paper; see that paper for details of the following assertions-specifically, see Theorem 2.4 and its proof. One can show 
that the controller-stopper problem has a continuation region given by {(x, y, z) : < x c / r (y, z) < x < xo(y, z)} for some 
functions < x c / r (y,z) < r/c < xo(y,z) with (y,z) £ R~. Thus, if x < x c / r (y,z), we have tfr(x,y,z) = (c/r)x, and if 
x > Xo(y,z), we have fi(x,y,z) = 1. Moreover, iff is non-decreasing and concave with respect to x on R + (increasing 
and strictly concave in the continuation region) and is the unique classical solution of the following free-boundary 
problem on [x c / r (y,z), Xo(y,z)]: 



$r(x c/r (y,z),y,z) = - r x c/r (y,z), ^ x (Xc/ r (y,z),y,z) = - r ; 
fi(x (y,z),y,z) = 1, &x(xo(y, z),y,z) = 0; 



in which 



-Cov 
-Civ 

X.2V 

M\v 
AI3V 



(m-y)v y + v Vyy, 
v V2n 



= -P12 



f(y,z) 



-A v — (r — A) x v x + 



1 



2\/(y,z) 



x 2 V, 



-P13 77 — r h{z)xv x: 

f(y,z) 

1 2 

g(z) v z + - h (z) v zz , 



= P23 V 



V2fc(z) 



Vyz, 
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(3.4) 

(3.5) 
(3.6) 

(3.7) 

(3.8) 

(3.9) 
(3.10) 



and 



NL 6 ' 6 = sup 

r 



+ r 



Pn 



(3.11) 



_ J(l -P 2 2 )(l -P 2 3 )-(P23 -Pl2Pl3) 2 
" (3> VtfAfe) -2 ; «Ax 



■P?2 



Because $ is concave with respect to x, we can express NIf' s as follows: 



NU* = -- v 2 (l - p 2 2 ) ^ - l -8h\z){\ - p 2 3 ) t^-v^J- h(z) (p 23 - P12P13) 



fyxyfyxz 



(3.12) 



5.2. Convex Legendre Dual of ij/ 

Since 0- is strictly concave with respect to x in its continuation region (which corresponds to wealth lying in 
[0, c/r]), we can define its convex dual *P by the Legendre transform: for (w,y, z) e D = {(w,y, z) e R 3 : w e [0, c/r]}, 



*F(w, y, z) = max (iA(jc, y, z) - wxj . 



(3.13) 



In this section, we show that the convex dual *P is the minimum probability of lifetime ruin; then, in the next section, 
we asymptotically expand ijt in powers of \[e and V^. 

Theorem 3.1. *P equals the minimum probability of lifetime ruin if/ on D, and the investment policy 71* given in 
feedback form by n* — n*(W*, Y t ,Z t ) is an optimal policy, in which W* is the optimally controlled wealth (that is, 
wealth controlled by n*) and the function n* is given by 



*, s H-r \p w 

n (w,y,z) = --z- — - P12 

f z (y,z) l/fym 



e f(y,z) i/fww 



<A..v ^ h(z) if/ K 

-P13 V<J 



f(y,z) i//» 



(3.14) 



in which the ri ght-ha nd side o/( |3.14| l is evaluated at (w, y, z). 

Proof. From (3.13 1, it follows that the critical value x* solves w = if/ x (x,y,z)', thus, given w, we have x* = I(w,y,z), 
in which / is the inverse function of \jr x with respect to x. Therefore, Y^y.z) = $r(I(w,y,z),y,z) - wl(w,y,z). By 
differentiating this expression of *F with respect to w, we obtain y,z) = ^/ x (I(w,y,z),y,z)I w (w,y,z) - I(w,y,z) - 
wl w (w,y,z) = -I(w,y,z); thus, x* = -*P vv (w,y,z). Similarly, we obtain (with w = {jf x (x,y,z)) the following expres- 
sions: 

1 



ij/ xx {x,y,z) 



fww(w,y,z) 
fy y {x,y,z) = *F v (w,y,z), 
fy t (x, y, z) = ^ziw, y, z), 
fiyy(x, y, z) = ^(w, y, z) (f/ xy (x, y, z) + *P, 7 (w, y, z), 
&zz(x, y> z) = ^wz(w, y, z) ^(jc, y, z) + ^(w, y, z), 
&ty(*, y, z) = ^ wy (w, y, z) iA.„(x, y, z), 
fefe y. z) = *P vl , z (w, y, z) iA.«(x, y, z), 



(3.15) 

(3.16) 
(3.17) 
(3.18) 
(3.19) 
(3.20) 
(3.21) 
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and 



fiy Z (x, y, z) = *P„3,(w, y, z) ^(w, y, z) <p xx {x, y, z) + ¥ yz (w, y, z). 



(3.22) 



By substituting x* = - ( P w (w,y,z) into the free-boundary problem for if}, namely (3.4i, one can show that *F uniquely 
solves the following boundary-value problem on D: 



ming Xfv{w, y, z) = 0; 
v(0,y,z) = 1, v(c/r,y,z) 



0. 



(3.23) 



where the differential operator ifi is given by 

Tfv = -A v + (rw + (// - r)P - c) v w + - (m - y) v y + 5g(z) v z 



1 



1 



1 



(3.24) 



+ ~ / (y, z)/3 v mv + - v z Vyy + -6 h L (z) v zz + p n f(y, z)/3v 
2 e 2 

+ Pnf(y,z)P ^fdh(z)v wz +p 2 3 V2v J^h(z)v yz . 
Observe that *F is strictly convex in w because if/ is strictly concave in x in its continuation region which corre- 



sponds to D in the original space. Since ¥ is strictly convex with respect to w, the optimal policy n* in ( 3.23 1 is given 



by the first-order necessary condition, which results in the expression in (3.14i. Now, using a standard verification 
theorem we deduce that *f is the minimum probability of lifetime ruin ifr. □ 



Theorem 



3.1 demonstrates the strong connection between if/ and if/, namely that they are dual via the Legendre 

z) 



transform. (As an aside, if we have i/r, we can obtain ip via Mx, y , z) = min H , (iff(w,y,z) + wx).) Therefore, if we have 
then we obtain the minimum probability of ruin iff via (3.13 ». More importantly, we get the optimal investment 



strategy n* via (3.14i. As a corollary to Theorem 3.1 we have the following expression for n* in terms of the dual 
variable x. 

Corollary 3.2. In terms of the dual variable x, tlie optimal investment strategy n* is given by n* — n*(X* , Y,,Z t ), in 
which X* is the optimally controlled process X, and 



K (x,y,Z) = ~-pr -X\p x 

f 2 (y,z) 



+ P12 A/ - 



e f(y,z) 



1 x <Jx h{z) f 

>Pxy +P13 77 T^x Z , 

/(y.z) 



(3.25) 



with the right-hand side of {3.25 I evaluated at (x,y,z). 



Proof. Let w = ij/ x (x, y, z) in ( 3. 14 1 and simplify the right-hand side via equations ( 3. 15 i-( 3.22 1 to obtain (3.25 1. □ 



4. Asymptotic Approximation of the Minimum Probability of Lifetime Ruin 

In this section, we asymptotically expand the Legendre transform of the minimum probability of ruin, in 
powers of yfi and y/S. (A parallel analysis of expanding the Legen dre transform of the value function of the utility 
maximization problem was carried out in |jonsson and Sircar| ( |2002| ).) We expand iff instead of \p because if one were 
to do the latter, then one would note that each term in the expansion solves a non-linear differential equation. The 
differential equation for the zeroth-order term has a closed-form solution; however, none of the differential equations 
for the higher-order terms does. What this fact implies is that to solve any of these non-linear differential equations, 
one would have to assume that it has a convex solution, determine the corresponding linear free-boundary problem for 
the concave dual, solve this free-boundary problem, then invert the solution numerically, as in equation ( |3. 1 3| >. Note 
that one would have to perform this procedure for each higher-order term in the expansion. 

By contrast, when we expand $r, each term solves a linear differential equation, as we show below. We explicitly 
solve these linear differential equations, then invert the approximation using p,13[ ) once to obtain an approximation 
for the minimum probability of lifetime ruin xfi. Note that the resulting approximation of if/ is not guaranteed to be a 
probability, that is, to lie in the interval [0, 1]; however, our numerical experiments show that this is not a problem for 
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the values of the parameters we consider. See Fouque et al. (2008) for an example of approximating a probability that 
solves a linear differential equation. 

To begin, expand if/ and the free boundaries in powers of V#: 



and 



(j/ = (£o + V<5 $1 + 5 xffz + ■ ■ ■ , 
x c / r (y,z) = x c/rfl (y,z) + VSx c / rt i(y,z) + Sx c/n2 (y,z) + ■■■ 

xo(y,z) = x ,o(y,z) + yfdx (u (y,z) + 6x Qt2 (y,z) + ■■■ . 



Insert the expression in (4.1 1 into Nlf' 5 in (4.12) to obtain the following expansion in powers of yd: 



NL 



1 2(t 2\ ^0,xy 
= --V 1 - Pn ) 



~ v { l -Pi2) h — 



2 ? \ 

fihxx ~ 2 lf/l,„ 

<P(),xx 



2 .... ,<Ao,*>.^0,. tZ 

vh(z)(p2i -pnPn) — ^ 



+ 0(6). 



Keeping terms up to Vd, we expand the free-boundary conditions in ( 3.4 1 as 

4io(x c ir&(y,z),y,z) + ^f5\x c i r ^(y,z)4/oJx cir fi(y,z),y,z) + 4>\ (x c/rfi (y, z), y, z)] 
= ~ (x c /r,o(y,z) + V# x c/rJ (y,z)) , 

$QA x cirfl(y,z),y,z) + ^6[x c / r ,i(y,z)iffo, X x(xc/ r fi(y,z),y,z) + &i,x( x c/r,o(y,z),y,z)] 



and 



<AoOo,oCy, z),y,z) + xo,i(y,z)^(x jo(y,z),y,z) + h(xo.o(y,z),y,z)] = 1, 
<fro,x(x ,o(y,z),y,z) + yf5[xo,i(y, z) fio,. xx (x (U) (y, z),y, z) + fii, x (xo,o(y,z),y,z)] = 0. 



(4.1) 
(4.2) 

(4.3) 



(4.4) 



(4.5) 

(4.6) 
(4.7) 
(4.8) 



We begin by approximating {j/Q and the free boundaries x c / r fl and xq q. Then , we use the boundaries x c / r ,o an d xo,o as 
fixed boundaries to determine ty\ . As one can see from equations (4.5 i-(4.8 1, this fixing of the boundaries introduces 
an 0( V(5)-error into in 0{ V<5)-neighborhoods of x c / r ,o and jto ; o. 

Terms of order 6° 



By inserting (4.1 i-(4.4i into (3.4i and collecting terms of order 6°, we obtain the following free-boundary problem: 

' cx + (I £ + ^ £i + - ^ (1 -p? 2 ) = 0; 

&o(x c/n0 (y,z),y,z) = - r x c/rfl (y,z), fto,x(x c / rt0 (y,z),y,z) = f, 
fro(x Qfi (y, z), y, z) = 1, $o,*(*o,oCy. z), z) = 0. 



(4.9) 
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Terms of order yd 

Similarly, by comparing terms of order Vtf and using x c / r fl and xo,o as fixed boundaries for tj/\, we obtain the 
following boundary-value problem: 



(i -Co + ^ -& + X 2 )<Ai + (Mi + ^ M 3 ) «Ao 

-:-(•-/-)((;: 

<Ai (* c / r ,o(y> z). y» z) = 0, 1A1 (^(xoCv, z), y, z) = 0. 



(4.10) 



Next, we expand the solutions of (4.9 1 and ( 4.10| > in powers of yfe: 

ij/o{x,y,z) = 4t Q fl(x,y,z) + yfe ^(x, y,z) + e^ 0i2 (x,y,z) + 

and 

iffi(x,y,z) = 4>]_fl(x,y,z) + VeiA u (x,y,z) + e^i t2 {x,y,z) + 
Similarly, expand the free boundaries x c / n o and xo$ in powers of -\fe: 

x c irfi{y,z) = x c/rfl fl(y,z) + yfex c/rA1 (y,z) + ex c/rA2 (y,z) + ■ 

and 

xo,o(y,z) = xo,o,o(y,z) + yfex 0A i(y,z) + exo,o,2(y,z) + ■■■ 



(4.11) 
(4.12) 

(4.13) 
(4.14) 



Substitute (4.1 1 1 and (4.12 1 into (4.9 1 and (4.10 1, respectively, and collect terms of the same order of yfe. As discussed 
earlier, we determine the free boundaries x c / rA o(y, z) and JCo,o,oCy> z) via a free-boundary problem for i^o.o; then, we use 
these boundaries as the fixed boundaries for i/^i and (^i q. 

Terms of order 1/e in \A.9) 



By matching terms of order 1/e in (4.9 1, we obtain the following: 



<Ao,o,-t 



or equivalently 



/ \ 1 It lit 2 \ v,u,xy „ 
(m - y) i// 0Ay + v iAo,0,rv- v V- ~ Pii) ^ = 0. 



<Ao,o,a- 



(4.15) 



(4.16) 



We, therefore, look for an $o.o independent of y; otherwise, 1^0,0 will experience exponential growth as y goes to ±oo 
( Fouque et al. 2000 , 2003 ). We also seek free boundaries Jc c / r ,o,o and xq,o,o independent of y. 



Terms of order 1 / yfe in (4.9 I 

By matching terms of order 1/ Ve in (4.9 1 and using the fact that (&o,o,y = 0, we obtain the following: 

Xo<Ao,i=0. (4.17) 

Therefore, we look for an $o,i independent of y; otherwise, 0-q i will experience exponential growth as y goes to +oo. 
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Terms of order e in (4.9 1 

By matching terms of order e in (4.9 1 and using the fact that i^o.Oj- = $a,\,y = 0, we obtain the following Poisson 
equation (in y) for 1^0,2: 

-Co *Ao,2 = -cx - £ 2 <Ao,o- (4. 18) 

The solvability condition for this equation requires that cx + £2 1^0,0 be centered with respect to the invariant distribu- 
tion N{m, v 2 ) of the Ornstein-Uhlenbeck process Y. Specifically, 



(cx + £2 iAo,o) = cx + (£ 2 ) iAo,o = 0, 
in which (■) denotes averaging with respect to the distribution N(m, v 2 ): 



<v> 



1 r° 

= VI 

\nv 2 J -00 



(y) e 2y 2 dy. 



(4.19) 



(4.20) 



In (4.43), the averaged operator (£2) is defined by 



1 / fi — r 
{£2) v = -A v - (r - A) x v x + - 



2 \cr,(z) 



x 2 v x 



(4.21) 



in which 0"«(z) is given by 



1 



crl(z) \f 2 (y,z) 



1 



Thus, we have the following free-boundary problem for i^o,o- 



(4.22) 



cx- A <Ao,o - (r - A) x {fr 0A x + s(z) x 2 tfr ,o,xx = ; 

&),oOc/r,0,o(z),Z) = 7 X c/rfl fl(z), tAo,0 ; x(* c /r,0,o(z),z) = 7 : 

&>,o(*o,o,o(z), z) = 1. ^o,o,^(^o,o,o(z)»z) = 0. 



(4.23) 



with s(z) = \ ■ The general solution of the differential equation in (4.23 i is given by 

foofoz) = D x (z)x B ^ + D 2 (z)x B ^ + -x, 



(4.24) 



in which 



and 



B 2 (z) 



1 



2s(z) 
1 



(r-A + s(z)) + J(r-A + s(z)) 2 + 



4As(z) 



tr- A+ v(:))- yir-.l + Mrir +4,!.s(:) 



> 1, 



<0. 



We determine D\ and D2 from the free-boundary conditions. 
The free-boundary conditions imply that 



(4.25) 



(4.26) 



Dl(z)X c / rfifi (z) Bi(Z> + D 2 (Z) X c / rfi ,o(z) B2(Z> + - Xc/r,0fi(z) = - Xcjrfifliz), 



D 1 (z)B I (z)x e/r , ,o(z) fil(zK1 +D2(z)B 2 (z)x c/rA o(z) 



B 2 fe)-1 



r r' 



Ddz)x QAl) (zf' (z) +D 2 (z)x 0A0 (zf 2{z> + -*o,o,o(z) = 1, 

r 



and 



Di(z) BiCz) *o,o,o(z) Bl(zM + 2 (z) B 2 (z) x 0A0 (z) 



B 2 (z)-l 



+ - = 0, 

r 



(4.27) 

(4.28) 
(4.29) 

(4.30) 



which gives us four equations to determine the four unknowns D\, D 2 , x c / r ,o,0< an d *0,o,o- Indeed, the solution to these 
equations is 

1 (c Bi(z)-1' 



Di(z) = - 



Bi(z)-l\r Bi(z) 
D 2 (z) = 0, 

X c /r,0.o(z) = 0, 



and 



It follows that 



B^z) r 



BAz)-l c 



<Ao,oO, z) = - 



1 /c Bi(z)-1 



Bi(z)-l\r B,( Z ) 



x. 



(4.31) 

(4.32) 
(4.33) 

(4.34) 
(4.35) 



Terms of order yfe in (4.9 1 

By matching terms of order -\[e in (4.9 1 and using the fact that 0o,O,y = tAo.i.v = 0, we obtain the following Poisson 
equation (in y) for 1^0,3 : 

-Co $"0,3 = -Xi «Ao,2 - X2 *Ao,i- 



As above, the solvability condition for this equation requires that 

(Xi ^0,2 + X2 <Ao,i ) = 0, 

in which 



It follows that $0,1 solves 



iAo,2(*, z) = X 1 (-« - X2 iAo,o) • 
(X2) (Ao,i = (XiXq 1 (cx + £ 2 iAo.o)) ■ 



Recall that we impose the (fixed) boundary conditions 1^0,1 C*r/r,o,o(z),z) = and 
fto,i(*o,o,o (z),z) = at x c/n0 ,o(z) = and x ,o,o(z) = ' c- 

From (4.38 1, it is straightforward to show that 1^0,2 can be expressed as follows: 

&Vi{x,y,z) = -£» 1 (z)B 1 (z)(B,(z) - l)x B ' (z) r,(y,z), 

in which 77 solves 

It follows that the right-hand side of (4.39 1 equals 



*(z). 



/(y.z) 



-p 12 Ou - r) v V2 Dife) U?Cz) (-Bi(z) " 1) ^ l(z) I 
= Pi 2 (/*-r) V2yD,(z)B?(z)(Bi(z)-l)^- (z) |%,z) 1 1- 



1 ( /j. - r 

7o^) 



-s(z) 



in which F is an antiderivative of l/f with respect to y; that is, 



F y (y,z) = 



1 



From (4.39 1 and (4.42 1, we obtain that i/Vo,i equals 

$o,i(x,z) = 5i(z)* Bl(z) + D 2 (z) x* 2 ® + A(z)x B ' (z) Inx, 
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(4.36) 

(4.37) 
(4.38) 
(4.39) 



(4.40) 



(4.41) 



(4.42) 



(4.43) 



(4.44) 



in which B\ and B2 are given in ( 4.25) > and ( |4.26 1, respectively, and A is given by 



A(z) = 



(2B l (z)-l)s(z)-(r-A) 



P(y,z) 



I ( fi — r 



2\f(y,z) 



-s(z) 



(4.45) 



The functions D\ and D2 are given by the (fixed) boundary conditions at x c / r ft,o(z) = and jco,o,o(z) = b^)-i ' c' fr° m 
which it follows that 



fro.i(x,z) = A(z)x Biiz) Jinx -In J 



«i(z) r 



= A(z)x B,(z) ln\x 



«i(z)-l c 
fli(z)-l £ 
Bi(z) ' r 



(4.46) 



Next, we focus on (4.10 1 to find #\o, after which we will approximate # by #0,0 + V^<Ao,i + V^*Ai,o- 
Terms of order l/e in ( |4.10[ ) 



By matching terms of order 1 fe in ( 4. 10 1, we obtain the following: 

-Co #1,0 =0, (4.47) 

from which it follows that \V\# is independent of y; otherwise, #10 will experience exponential growth as y goes to 
+00 dFouque et aL][2003> . 



Terms of order 1/ yfe in (4.10l 

By matching terms of order 1 / -\fe in ( 4. 10 1 and using the fact that #1,0,3) = 0, we obtain the following: 

-Co #u=0. (4.48) 

Therefore, we look for an #11 independent of y\ otherwise, \ will experience exponential growth as y goes to +00. 



Terms of order e° in (4.10 1 

By matching terms of order e° in (4.10 1 and using the fact that #1,0,31 = #U,y = 0, we obtain the following Poisson 
equation (in y) for \jf\£. 



fjL — r 

-Co #1,2 = --£2 #1,0 + P13 -t. — r Hz) x if/o,o,xz- 

/(y.z) 



The solvability condition for this equation requires that 



-£2 #1,0 + P13 -y, — c Kz) x #o,o,.vz ) = 0, 

/(y.z) 



or equivalently, 



(£2) #1,0 =Pi 



/(y,z) 



/i(z)x0- o> o,«- 



(4.49) 



(4.50) 



(4.51) 



with boundary conditions o(jc c /, oo(z),z) = and #1 o( x ooo(z),z) = at the boundaries x c / r00 (z) = and x oo(z) 

Bi(z)-1 ' c- 



Bl(z) ■ '- . It follows that 0-1,0 is given by 



ho(z)=x Blb) ln|x- 



in which A! and A 2 are 



Ai(z) = 



#l(z)-l £ 



ffi(z) 



Ai(z)+A 2 (z)lnbc 



«i(z) r 



Bi(z)-1 c 



»2(z)i(z) 



(2Bi(z) - 1) s(z) -(r-A) [(2Bi(z) - 1) s(z) - (r - /I)] 2 ' 
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(4.52) 



(4.53) 



and 



A 2 (z) = 



1 



H 2 (z) 



2 (2Bi(z)-l)jCz)-(r-A)' 



with H\ and H 2 functions of z defined by 

H l (z) + H 2 (z) lnx 
(i-r 



-pnh(z) 



B'Az) /Bi(z)-1 c 



/(y,z)/Bi(z)-l\ Bi(z) 



B,(z) 



1 +Bi(z) In be 



(4.54) 



(4.55) 



B!(z)-l\r Bi(z) 



B](z) 



4.7. 77ie Approximation of the Probability of Lifetime Ruin and the Optimal Investment Strategy 
Combining (4.35 i, (4.46 1, and (4.52 1, we obtain the following approximation of if/ 

iff e ' 6 {x,z) = ifrofl(x,z) + yfetfr 0A (x,z) + ^6ij/ U ){x,z) 
1 (c Bi(z)-1 



Bxiz) 



y/lA(z)x Bl(z) ln jc- 



VSjc Blfe) InU- 



gi(z)- 1 £ 
Bi(z) ' r 
Bi(z)-l c_ 
Bi(z) ' r 



(4.56) 



A 1 (z)+A 2 (z)mx- 



Bi(z)-l c 



in which A, A u and A 2 , are specified in (4.45 i, (4.53 i, and (4.54 1, respectively. 

We also a pprox imate the dual of the optimal investment strategy up to the first powers of yfe and Vc5, as we did 
for {j/. Using (3.25 i, we obtain 



ft £ ' e (x,z) 



f 2 (y,z) 
- W- 



pi — r v V2 ~ 

■ X l/r , 1 ,.xx + P 1 2 -J7— ~ ^0,2,) 



fHy,z) 



f(y,z) 



p-r 



f(y,z) 



xfy\,a,xx+ P\3 



h(z) 

f(y,z) 



(4.57) 



Given w e IR+, we solve for x using w — if/x d (x,z)- Then, we let t]/ €,s (w, z) := (j/ € ' 6 (x, z) - xw, thereby performing 
the calculation in equation ( |3. 1 3) >. We also denote by n e - d the function that satisfies n e,5 (yv,z) :— Tt e,s (x,z). Note that 
the resulting approximation of tfr is not guaranteed to be a probability; however, this is not a problem in the numerical 
experiments we consider in the next section. 



5. Numerical Solution using the Markov Chain Approximation Method 

In this section, we describe how to construct a numerical algorithm for the original optimal control problem 
directly using the Markov Chain Approximation Method (MCAM); see e.g. (Kushne r and Dupuisl^OOlUKus hner, 
1999) . For the ease of presentation, we will describe the numerical algorithm only when the fast scale volatility factor 
is present. In what follows p will denote the correlation between the Brownian motion driving the stock and the one 
driving the fast factor, that is, p = p\ 2 . 

Let us fix an /i-grid, that is, a rectangular compact domain G h c R 2 with the same spacing h in both directions. 
We choose an initial guess (on this grid) for a candidate optimal strategy. Denote this strategy by ji. Then, our goal is 
to create a discrete-time Markov chain (^')«>o that lives on G h and that satisfies the local consistency condition 

Eg[A&,] = b{x,n)Af\x,n) + o(Af /! ), 
Cov'%[A{ h n+1 ]=A(x,7r)Af\x,7T) + o(At h ), ' 

in which Af„ + i = £„ + i - £„, and b and A denote the drift and the covariance of the vector X, = {W,, Y,), respectively. 
(The Markov chain is constructed to approximate this vector in a certain sense.) E'jJ denotes the expectation, given 
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that the state of the Markov chain at time n is x. In \5.\\ the quantity At h (called the interpolation interval) is to be 
chosen so that it goes to zero as h — > 0. We also do not want this quantity to depend on the state variables or the 
control variable. 

Since G h is a compact domain, we impose reflecting boundary conditions at its edges. (Natural boundaries exist 
for W(t), specifically and the safe level e . However, Y, lives on an infinite region.) For example, we choose the 
transition probabilities to be p"' ((w,y), (w,y - h)) = 1, when y is as large it can be in G h and for all we [0,f]. 

5.1. Constructing the Approximating Markov Chain 



5.1.1. Whenp = 0. 

Denote a - j,/? = v y|. We obtain the transition probabilities of the Markov chain % h as 

/3 2 /2 + ha(m -y)± 



P K - n ((w,y),(w,y + h)) 



p n ' h ((w,y),(w±h,y)) = 
y),(w,y)) 



and choose the interpolation interval to be 



Qh 

(f(y)n(w,y)) 2 l2 + hip, - r^wj) 1 + h(rw - cf 



Q h -Q"- h (w,y) 



Af 



Q" 



in which 



and 



Q Kj '(w,y) = (nfiy)) 1 + j3 l + h\a(m - y)\ + h\(p - r)n(w,y)\ + h\rw - c\, 



Q h = max &*(w,y), 

(w,y,7t) 



in order to satisfy the local consistency condition. Here a ± = max{0, ±a) 



(5.2) 



5.1.2. WhenptO. 

In this case a convenient transition probability matrix solving the local consistency condition is 



p x - h ((w,y),(w,y±h)) 



p"' h ((w,y),(w±h,y)) = 



(1 - p z )p z /2 - \pn(w,y)\l3f(y)/2 + ha(m - yf 

(f(y)n(w,y)) 2 /2 - \pn(w,y)\p-f(y)/2 + h{p- r)n{w,yT + h(rw - c)± 



p x - h ((w,y), (w + h,y + hj) = p Kjl ((w,y), (w - h,y - h)) 



p K ' h ((w,y), (w + h,y - Kj) = p"' h ((w,y), (w - h,y + h)) 
p"- n ((w,y),(w,y)) = , 



Qh 

(pn( W ,y)T/3f{y) 
2Q h 

(pn( W ,y))-l3f(y) 
2Q h 



where 



Q nJl (w,y) = (nf(y)f + f - \pn( W ,y)\f3f(y) + h\a(m - y)| + h\Ql - r)n(w,y)\ + h\rw - c\. 



(5.3) 



For values of \p\ close to 1, the transition probabilities may be negative. The positiveness of these probabilities is 
equivalent to the diagonal dominance of the covariance matrix A = (0,7). (Recall that we call A diagonally dominant 
if an — 'Eijjti \ a ij\ > 0, V/.) The construction of an approximating Markov chain when some of the expressions in 
(|5.3|> are negative will be discussed next. 
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5.1.3. When p — 1 and some of the transition probabilities in ( |5.3| l are negative. 

We accomplish the construction of the approximating Markov chain in two steps, following Kushner ( 1999): 
(i) Decomposition. As in |Kushner and D upuis (2001 ) Sections 5.3 and 5.4, we decompose X into separate components 
and build approximating Markov chains to match each component. Then, we combine the transition probabilities 
appropriately to obtain the approximating Markov chain for X itself. 

Let X = X m + X <2 \ in which 



dX, 



dX) 



(2) 



"Tr- 

rW, - c + (fi- r)n, 
a(m — Y,) 



dt. 



(5.4) 
(5.5) 



Since p = 1, we take B l = B 2 . Suppose that the form of the locally consistent (with dynamics of and Z (2) , 
respectively) transition probabilities and interpolation intervals are 



P{ (x,x) = 



n.h/ -\ 

n x (x,x) 



Q 



1 Q'l 



\ ( 



- Q h 2 



for some n n ^ h (x, x), n^ h {x, x), and appropriate normalizers Q h v Q 1 ^. Then, the following transition probabilities and the 
interpolation interval are locally consistent with the dynamics of X 



p"' (x,x) 



n\' h (x,x) + hn^ h {x,x) 



Af' = ^z- 



(5.6) 



Since it is easier, we first provide the expression for p 2 



P{ ((w,y),(w,y + h)\n) = 



P2 h ((w,y),{w + h,y)\n) = 



pl' h ((yv,y),(w,y)) 



(p. - r)7r(w,y) ± + (rw - c) 1 

z ^ 

Q2-Q n 2 '\w,y) 



(5.7) 



where 



Q n ? '\w,y) = a\m -y\ + {p- r)\n(w,y)\ + \rw - c\. 



The computation of p"' h is more involved. This is the subject of the next step. 



(ii) Variance control. System (5.4 1 is fully degenerate; that is, the corresponding co variance matrix A is not diagonally 
dominant. The previous technique for building a Markov chain does not work. Instead, we will build an approximating 
Markov chain by allowing the local consistency condition to be violated by a small margin of error. 

If (<xi , 1x2) = (qk\ , qkj) for some constant q and integers k\ , k 2 , we could let the transition probability to be p h (x, x± 
(hk\,hk2)) =1/2 and the interpolation interval to be At h = h 2 /q 2 , and we would obtain a locally consistent Markov 
chain. This is not possible in general. For an arbitrary vector (<t\ , crj), we can find a pair of integers k\ (x, n), k 2 (x, ri), 
and a real number y(x, n) e [0, 1], such that 



0"i(x,7r) 
0-2 



q(x, 7T) 



k\{x,n) 
k 2 (x, 7t) + y(x, n) 
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Since the Markov chain is constrained to the grid G h , we can only approximately let it move in the direction of 
(ctj , o- 2 ) T . We choose 

p n > h (x,x±h(k 1 ,k 2 ) T )=pi/2, 

p n \x,x±h{k u k 2 + l) r ) = p 2 /2, D ' 8) 

in which p\ + p% = 1, and /?i and /?2 will be appropriately chosen in what follows. The mean and the covariance of 
the approximating chain is 

E^[A£W)] = 0, 



E*J[A£*(jc, n)A?(x, n) T ] = h 2 C(x, n), ' ? ' 9 ' 



where 



,., , , k\ k x k 2 \^ I k\ h(k 2 + l) 
C(x,n)=pA, ) , 2 \ + p2 l 



k x k 2 k\ ) ^{hika + l) (k 2 + l) 



k\ k x (k 2 +p 2 ) 
(k 2 + p 2 ) k\ + 2pk 2 + p 

We choose the interpolation interval to be Af' h (x, n) = h 2 /q 2 . On the other hand 



(5.10) 



1 \ m \i 2 / k \ hiki+y) 
a(x,7t) = A(x,n) q = , ,, , . 2 
y \k\{k 2 +y) {k 2 +yf 

and we see that if we pick p 2 = y, then C\\ = ci\\ and C\ 2 = a\ 2 match, but we violate the local consistency condition 
by 

( 1 ^ 



c 22 ~ a 22 _ y(l -y) _ () 
a 2 (k 2 + y) 2 



(5.11) 

We will choose k 2 sufficiently large so that the local consistency condition is almost satisfied, and the numerical noise 
in ( |5.1 l| l is significantly reduced. 



5.7.4. The case when p e (— 1, 1) and some of the transition probabilities in \53) are negative 
We will decompose the state variable into three components: 

^_lW\_lrW,-c + {p-r)n,\ Jn,f{Y,)\ l jQ WdfiA 
dX '-\dY t j-\ a(m-Y t ) J dt + [ P p j dB < + [o p ^1 l^j ' ( } 

that is, a drift component, a fully degenerate noise component, and a noise component with diagonally dominated 
covariance matrix. We can build an approximating Markov chain for each component separately and then combine 
them as discussed above. 

5.2. Approximating the Probability of Ruin and Updating the Strategy 
We solve the system of linear equations 

V*>\x) = e-"** Yj P n ' h(x > x) v "' h (x)> (5- 13) 

xeC 1 ' 

with boundary conditions V K,ll (Q,y) = 1 and V K ' h (c/r,y) = 0. This is the dynamic programming equation for a 
probability of ruin problem when the underlying state variable is the Markov chain £*. In the next step, we update our 
candidate for the optimal strategy. For convenience, denote V' 1 ' by V. In the interior points of the grid 

, Hp. - r)[V(w + h,y)- V(w,y)] + (J3/2)pf(y) [V(w + h,y + h) + V(w,y - h) - V(w + h,y - h) - V(w,y + h)] 

7T( W. V ) = . 

f 2 00 [ V(w + h,y)+ V(w -h,y) -2 V(w,y)] 
On the wealth dimension boundaries of the grid, we let n{c/r,y) = and 

hQx - r)[V(h,y) - V(0,y)] + (fi/2)pf(y) [V(h,y + h) + V(0,y - h) - V(h,y- h) - V(0,y + h)] 
f 2 (y)[2V(0,y) - 5V(h,y) + 4V(2h,y) - V(3h,y)] 
The updates of the optimal strategy for the maximum and minimum values of y are similar. 
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Iteration 

Once the optimal strategy is updated, we go back and update the transition probabilities and solve the system of 
linear equations in ( |5.13[ ) to update the value function. This iteration continues until the improvement in the value 
function is smaller than an exogenously picked threshold. 

Two Technical Issues 

• The initial guess of the optimal strategy is important. For p = 0, we take the initial strategy as the one in 



constant volatility case, where the closed-form solution is available in Young (2004). Forp + 0, we take the 
final strategy computed from zero-correlation case (p = 0) as the initial guess. This initial guess makes the 
algorithm converge fast. 

• For p + 0, the covariance matrix of the wealth process and volatility factor, in general, does not satisfy the 
diagonal dominance condition. The problem is more serious for the slow factor, since its variance is of the 
order of 6, and the numerical noise using '"variance control'" is far greater. To solve this issue we perform a 
"scale adjustment" to increase the variance of the factor. For example, if we define Z f = 100Z,, then the dynamic 
of the system becomes 

dZ, = 6(1 00m - Z,)dt + 100 Vtf V2v 2 dBf ) . 

when g — (m - z) and h = V2v. The new system is mathematically equivalent to the original one, but with a 
much bigger variance; thus, the numerical noise in variance control is much smaller. Note that the scheme here 
is equivalent to choosing a different grid sizes for the volatility and wealth dimensions. 

6. Numerical Experiments 

In order to conduct our numerical experiments we will take the dynamics of the slow factor in ( j2.7| i to be 

dZ, = 8{m - Z,)dt + Vtf V2vdB| 3) , Z = z. 



We let f(y, z) = exp(— y) or f(y, z) = exp(-z) in ( |2.2| l, depending on whether we want to account for the fast volatility 
factor or the slow volatility factor in our modeling. We will call 1 je or 5 the speed of mean reversion. We will take 
the correlations between the Brownian motions driving the volatility factors and the stock price to be p = pn = pn- 
The following parameters are fixed throughout this section: 

• r — 0.02; the risk-free interest rate is 2% over inflation. 

• /i = 0.1; the expected return of risky asset is 10% over inflation. 

• c = 0. 1 ; the individual consumes at a constant rate of 0. 1 unit of wealth per year. 

• A — 0.04; the hazard rate (force of mortality) is constant such that the expected future lifetime is always 25 
years. 



1.364 and v = 0.15, so that the harmonic average volatility, which we will denote by cr m = -\/l/E[l// 2 (T)] 



Vl/E[e 2y ] = e-™-" 2 = 0.25, in which Y is a normal random variable with mean m and variance v . The distri- 
bution of this random variable is the stationary distribution of the process (T/) f >o; see ( |4.22| l. [Note that o~ m is 
very close in value to E[/(F)] = E[e~ F ] = e -" 1+v2 / 2 = 0.26.] 

In our numerical procedure we use a bounded region for Y and impose reflecting boundary conditions. However, 
f(Y t ) is not bounded and not bounded away from zero. On the other hand, the invariant distribution of the process Y 
is normal with mean 1.364, and variance 0.15 2 . So it is with very small probability that Y, is negative or very large. 
Therefore, the fact that f{Y t ) is not bounded or bounded away from zero does not affect the accuracy of our numerical 
work in a significant way. 
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Observation 1 

We give a three-dimensional graph of the minimum probability of ruin and the optimal investment strategy in 
Figure [T] which are computed using MCAM. Here the speed of mean reversion is 0.5, p — 0, and only one factor 
is used. In our experiments we observed that the optimal strategy n* is positive (no-shortselling). As expected we 
observe that w — > i[/(w,y) is convex and decreasing. Note that f(y) — > if/(w,y) is increasing. Also, f(y) — > n*(w,y) 
is decreasing; however, it is not necessarily true that w — > n*(w,y) is decreasing. The latter behavior depends on the 
value of y. 

The probability of ruin does not depend on the sign of the correlation, p, between the Brownian motions driving 
the stock and the one driving the volatility. The larger the magnitude of p, the larger the probability of ruin. However, 
the minimum probability of ruin is quite insensitive to the changes in p; see Figure [2] 

Observation 2 

We compare the optimal investment strategy n*{w,y,z) in ( |3.14| i to 

fx- r c - rw 

n(w; cr) = — — — , 

cr 1 (p - \)r 

in which 

p = — [(r + A + s) + y/{r + A + s) 2 - 4rA\ , 

and 

When we want to emphasize the dependence on cr, we will refer to p as p{cr). Young (2004) showed that the strategy 
jt is optimal when the volatility is fixed to be cr. 

If only the fast factor is present and the speed of mean reversion is 250 (e = 0.004), then tf/ e,s in ( |4.56| l can be 
expressed as 

iA £ 0) = 0o,oO) 

whose inverse Legendre transform is 

<A>;cr m ) = (l--w) , (6.1) 

which is exactly the minimal probability of ruin if the volatility were fixed at cr m . Therefore, it is not surprising that 
for very small values of e, the minimum probability of ruin i//(w,y), calculated using MCAM, can be approximated by 
( |6.1| i; see Figure [3]-a. In our numerical calculations and in ( |6.1| l, we observe that the minimum probability of ruin ifr 
does not depend on its second variable. This result is intuitive, since when only the fast factor is present whatever the 
initial value of cr is, the volatility quickly approaches its equilibrium distribution (which is normal with mean cr m ). 
In fact 7To(w; o~ m ) practically coincides with the optimal investment strategy it* , which is computed using MCAM; see 
Figure [3}b. 

The most important conclusion from Figure|3]-b is that it is not necessarily true that the optimal investment strategy 
when there is stochastic volatility is more or less than the optimal investment strategy when the volatility is constant. 
Comparing n(w; cr) and n*(w, - ln(cr)) for different values of cr, we see that tt*(w, - ln(cr)) < n(w; cr) for larger values 
of cr, whereas the opposite inequality holds for smaller values of cr. The investment amount decreases significantly as 
the volatility increases. 

If only the slow factor is present and the speed of mean reversion is 0.02, then 



ip s (w,z) = \l- -w\ , (6.2) 



\p(e- z ) 

V 

c 

approximates the minimum probability of ruin if/(w,z), which we calculate using MCAM, quite well; compare 
ip(w, — ln( cr)) and if/ s (w, -ln(cr)) for different values of cr in Figure |4[a. We also compare n(w; cr) and n*(w, -ln( cr)) 
for several values of cr and draw the same conclusions as before. Also note that the optimal investment strategy is not 
necessarily a decreasing function of wealth. 

When we take the speed of mean reversion to be 0.2 (medium speed), then the probability of ruin starts diverting 
from what j6.l\ or ( |6.2| i describes; see Figure |5}a. As to the comparison of the optimal investment strategy with 
7r(w; cr), the same conclusions can be drawn; see Figure [5]-b. 
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Observation 3 

We compare the performance of several investment strategies in the stochastic volatility environment. Let <xo be 
the initial volatility. We denote by n M the strategy when one only invests in the money market. The corresponding 
probability of ruin can be explicitly computed as if/ M (w) = (1 - c/rw) We will also denote 7f"(w) = Tf(w;o~o), 

7T b = 7r(w, cr m ), and 

u — r c — rw 

*(>W) = 7^7 7T- 6 - 3 ) 

f z (y,z) (p- y)r 



Let n e (w) denote the approximation to the optimal strategy we obtained in Section 4.1 when we only use the e- 
perturbation. Similarly, let ^(w, z) be the approximation to the optimal strategy when we only use the ^-perturbation. 

We obtain the probability of ruin corresponding to a given strategy n by solving the linear partial differential 
equation Tfv = 0, (see ( |3.24| > for the definition of the differential operator D") with boundary conditions v(0, y, z) — 
and v(c/r,y,z) = 1. (This computation uses the MCAM without iterating.) 

In Figure [6]|8] we observe that the performance of if and if are almost as good as the optimal strategy n*. (Here 
we are considering a medium mean reversion speed. When the mean reversion speed is much smaller, then n s would 
be a better investment strategy.) Moreover, their performances are robust, in that, they do not depend on the initial 
volatility cr . This should be contrasted to if and tt. The former performs relatively well when cr is small, whereas 
the latter performs better when <tq is large. When <tq - cr m , all strategies perform as well as the optimal strategy. 
Also, observe that for wealthy or very poor individuals the choice of the strategy does not matter as long as they invest 
in the stock market. The difference is for the individuals who lie in between. 

As a result, we conclude that if the individual wants to minimize her probability of ruin in a stochastic volatility 
environment, she can still use the investment that is optimal for the constant volatility environment. She simply needs 
to update the volatility in that formula whenever the volatility changes significantly. 



References 

Bayraktar, E., Moore, K.S., Young, V.R., 2008. Minimizing the probability of lifetime ruin under random consumption. North American Actuarial 
Journal 12(4), 384-400. 

Bayraktar, E., Young, V.R., 2007a. Correspondence between lifetime minimum wealth and utility of consumption. Finance and Stochastics 11, 
213-236. 

Bayraktar, E., Young, V.R., 2007b. Minimizing the probability of lifetime ruin under borrowing constraints. Insurance: Mathematics and Economics 
41, 196-221. 

Bayraktar, E., Young, V.R., 2008. Minimizing the probability of rain when consumption is ratcheted. North American Actuarial Journal 12(4), 
428^142. 

Bayraktar, E., Young, V.R., 2009. Optimal deferred life annuities to minimize the probability of lifetime ruin. North American Actuarial Journal 
13(1), 141-154. 

Bayraktar, E., Young, V.R., 201 1. Proving regularity of the minimal probability of rain via a game of stopping and control. Finance and Stochastics 

(in press) Available at http://arxiv.org/abs/0704.2244. 
Fouque, J.R, Papanicolaou, G., Sircar, K.R., 2000. Derivatives in Financial Markets with Stochastic Volatility. Cambridge University Press, 

Cambridge. 

Fouque, J.R, Papanicolaou, G., Sircar, K.R., Solna, K., 2003. Multiscale stochastic volatility asymptotics. SIAM J. Multiscale Modeling and 
Simulation 2, 22^12. 

Fouque, J.R, Wignall, B., Zhou, X., 2008. First passage model under stochastic volatility. Journal of Computational Finance 11, 43-78. 
Jonsson, M., Sircar, K.R., 2002. Partial hedging in a stochastic volatility environment. Mathematical Finance 12, 375^109. 

Kushner, H.J., 1999. Consistency issues for numerical methods for variance control, with applications to optimization in finance. IEEE Transactions 

on Automatic Control 44, 2283-2292. 
Kushner, H.J., Dupuis, P., 2001. Numerical methods for stochastic control problems in continuous time, volume 24. Springer Verlag, New York. 

second edition. Stochastic Modelling and Applied Probability. 
Milevsky, M.A., Moore, K.S., Young, V.R., 2006. Asset allocation and annuity-purchase strategies to minimize the probability of financial rain. 

Mathematical Finance 16, 647-671. 

Moore, K.S., Young, V.R., 2006. Optimal and simple, nearly optimal rules for minimizing the probability of financial ruin in retirement. North 

American Actuarial Journal 10(4), 145-161. 
Young, V.R., 2004. Optimal investment strategy to minimize the probability of lifetime rain. North American Actuarial Journal 8(4), 105-126. 



18 



0.9 
0.8 




(a) Minimum probability of ruin 
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Figure 1: Minimum probability of ruin and optimal strategy computed by MCAM. Speed of mean reversion= 0.5. 
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(a) Probability of ruin at <tq=0.6 
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(b) Probability of ruin at <tq= (T flJ =0.25 
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(c) Probability of ruin at <to=0.1 

Figure 2: Variations of the minimum probability of ruin with respect to p. Speed of mean reversion^ 0.5. 



20 




wealth w 



(a) Minimum probability of rain 




wealth w 
(b) Optimal investment strategy 

Figure 3: Stochastic volatility versus constant volatility environment. Speed of mean reversion=250. 
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(a) Minimum probability of ruin 




(b) Optimal investment strategy 

Figure 4: Stochastic volatility versus constant volatility environment. Speed of mean reversion=0.02. 
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(a) Minimum probability of ruin 




(b) Optimal investment strategy 

Figure 5: Stochastic volatility versus constant volatility environment. Speed of mean reversion 
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wealth w 

Figure 6: Performance of the investment strategies described in Observation 3 (of Section 6). Speed of mean 
reversion=0.2. Correlation p = 0.5. Initial volatility cr () = 0.6. 




wealth w 

Figure 7: Performance of the investment strategies described in Observation 3 (of Section 6). Speed of mean 
reversion=0.2. Correlation p - 0.5. Initial volatility cr () = <x m = 0.25. 
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Figure 8: Performance of the investment strategies described in Observation 3 (of Section 6). Speed of mean 
reversion=0.2. Correlation p = 0.5. Initial volatility o- =0.1. 
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